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Abstract 

The quantization of Lorentzian or Euclidean 2+1 gravity by canonical methods 
is a well-studied problem. However, the constraints of 2+1 gravity are those of a 
topological field theory and therefore resemble very little those of the corresponding 
Lorentzian 3+1 constraints. 

In this paper we canonically quantize Euclidean 2+1 gravity for arbitrary genus 
of the spacelike hypersurface with new, classically equivalent constraints that maxi- 
mally probe the Lorentzian 3+1 situation. We choose the signature to be Euclidean 
because this implies that the gauge group is, as in the 3+1 case, SU(2) rather 
than SU(1, 1). We employ, and carry out to full completion, the new quantization 
method introduced in preceding papers of this series which resulted in a finite 3+1 
Lorentzian quantum field theory for gravity. 

The space of solutions to all constraints turns out to be much larger than the 
one as obtained by traditional approaches, however, it is fully included. Thus, by 
suitable restriction of the solution space, we can recover all former results which 
gives confidence in the new quantization methods. The meaning of the remaining 
"spurious solutions" is discussed. 



1 Introduction 

The canonical quantization of 2+1 (pure) gravity is a well studied problem and the lit- 
erature on this subject is extremely rich (see ]I| and references therein). It may appear 
therefore awkward to write yet another paper on this subject. 

The point of this paper is to quantize 2+1 gravity by starting with a new Hamiltonian 
(constraint) rather than the one that imposes flatness of the connection (see, for instance, 
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||). Therefore, we are actually dealing with a new field theory. The reason why we still 
can call this theory 2+1 gravity (although in the Euclidean signature) is because classically 
both theories are equivalent, it is in the quantum theory only where discrepancies arise. 

The motivation to study this model comes from 3+1 gravity : In a new method 
is introduced to quantize the Wheeler-DeWitt constraint for 3+1 Lorentzian gravity and 
one arrives at a finite quantum field theory. It is therefore of interest to check whether 
that quantum theory describes a physically interesting phase of the full theory of quantum 
gravity. One way to do that is to apply the formalism to a model system which maximally 
tests the 3+1 theory while being completely solvable. 

It is often said that 2+1 gravity in its usual treatment as for instance in |2| is such a 
model which tests the 3+1 theory in various technical and conceptual ways. The author 
disagrees with such statements for a simple reason : 

The constraints of usual 2+1 gravity and of 3+1 gravity are not even algebraically similar. 
Thus, one has to expect that the resulting quantum theories are mutually singular in a 
certain sense. We will find that this expectation turns out to be correct. 
One can partially fix this by studying Euclidean 2+1 gravity to test Lorentzian 3+1 gravity 
because then the two gauge groups (577(2)) coincide, in the Lorentzian signature the gauge 
group of 2+1 gravity would be 577(1,1). However, this is not enough : while now the 
Gauss constraints of both theories generate the same motions the rest of the constraints 
are still very different with respect to each other. More precisely, the 2+1 remaining 
constraint says that the connection A is flat, that is, its curvature F vanishes. Thus it 
does not involve the momenta E conjugate to the connection at all. The situation in the 
3+1 theory is very different : here we have as the remaining constraints a constraint that 
generates diffeomorphisms and the famous Wheeler-DeWitt constraint. Both constraints 
depend on the momenta, the Wheeler-DeWitt constraint even non-analytically. In M the 
authors propose to quantize the constraints FE = FEE = 0. However, this has never 
been done in the literature, one reason being that the FEE constraint is as difficult to 
quantize as in the the 3+1 case. Moreover, the two constraints FE = FEE = are 
equivalent to the F = constraint only when the two-metric q is non-singular, that is, 
det(g) > and therefore it is no surprise that the two theories are not even classically 
equivalent as was shown in || (for the theory defined by the F = constraint the condition 
det(g) > is put in by hand in order to have Euclidean signature). 

In this paper we are using the constraints FE = FEE/ y / det(g) = 0. There are several 
reasons that speak for this choice : 

First of all these constraints are at least classically completely equivalent to the F = 
constraints because clearly they make sense only when det(g) > 0. In fact we will show 
that there is a field dependent non-singular map between the Lagrange multipliers of the 
two theories which map the two sets of constraints into each other. 
Secondly, they are just as in the 3+1 theory non-analytic in E (because det(g) is a func- 
tion of E) and so will test this feature of the 3+1 theory as well. In particular, both 
constraints are densities of weight one and only constraints of this type have a chance to 
result in densely-defined diffeomorphism covariant operators as argued in ||. 
Thirdly, these constraints are maximally in analogy to all the 3+1 constraints. 

The plan of the present paper is as follows : 

In section 2 we review the classical theory of Euclidean 2+1 gravity and outline our 
main strategy of how to arrive at a well-defined Hamiltonian constraint operator. 
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In section 3 we review the necessary background information on the the mathemati- 
cal tools that have been developed for diffeomorphism invariant theories of connections. 
Those Hilbert space techniques are identical for the 2+1 and 3+1 theory so that we have 
one more reason to say that the model under consideration tests the 3+1 situation. Also 
we need to construct a volume operator which as in the 3+1 theory plays a key role in 
the regularization of the (analog of the) Wheeler-DeWitt constraint operator. The 2+1 
volume operator turns out to be much less singular than the 3+1 operator which has some 
important impact on the regularization of the constraint operators. 

In section 4 we regularize the Wheeler-DeWitt operator. Many of the details are 
exactly as in the 3+1 theory although there are some crucial differences coming from the 
lower dimensionality of spacetime and also from the different singularity structure of the 
volume operator. 

In section 5 we perform various consistency checks on the 2+1 Wheeler-DeWitt oper- 
ator obtained, in particular whether it is a linear, covariant and anomaly-free operator. 

In section 6 we construct the full set of solutions to all constraints. It his where 
we encounter, besides reassuring results that give faith in the programme started in M, 
several surprises : 

• The quantum theory admits solutions which correspond to degenerate metrics. This 
happens although classically such solutions do not exist given our constraints. This 
should not be confused with the situation in || because there degenerate metrics 
are allowed even at the classical level. 

• We find an uncountable number of rigorous distributional solutions to all constraints 
which reveal an uncountable number of quantum degrees of freedom just as in any 
field theory with local degrees of freedom. This is in complete contrast to the usual 
treatment via the F = constraints which results in a topological quantum field 
theory with only a finite number of degrees of freedom. 

• The space of solutions contains the solutions to the quantum F = constraints 
as a tiny subspace. This subspace of solutions can be equipped with an inner 
product which is precisely the one that one obtains in traditional approaches. This 
is reassuring that our methods lead to well-established results and do not describe 
some unphysical phase of the theory. 

• The huge rest of the solutions cannot be equipped with the inner product appropriate 
for the F = constraints because they do not correspond to measurable functions 
with respect to the corresponding measure. However, there is another natural inner 
product available with respect to which they are normalizable. This inner product 
is likely to be the one that is appropriate also for the physically interesting solutions 
of the 3+1 constraints. The solutions to the F = constraint in turn are not 
normalizable with respect to this second inner product. Thus as expected, the two 
sets of constraints have solution spaces which lie in the same space of distributions 
but they cannot be given the same Hilbert space topology. It is in this sense that 
the quantum theories are mutually singular. 

In section 7 we conclude with some speculations of what the present paper teaches us 
for the 3+1 theory with regard to the solutions that are spurious from the point of view 
of the F — constraint. 
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In the appendix we compute the spectrum of the 2+1 volume operator for the simplest 
states. 

Throughout the paper we mean by the wording "2+1 or two-dimensional" always 2+1 
Euclidean gravity while by "3+1 or three-dimensional" we always mean 3+1 Lorentzian 
gravity. 

2 Classical Theory 

Let us start by reviewing the notation (see, for instance, 0). 

We assume that the three-dimensional spacetime is of the form M = RxE where £ is a 
two-dimensional manifold of arbitrary topology, for instance, a compact, connected two- 
dimensional smooth manifold, that is, a Riemann surface of genus g or an asymptotically 
flat manifold. Let e l a be the co-dyad on £ where a,b,c, .. = 1, 2 denote tensor indices and 
i,j,k,.. = 1,2,3 denote su{2) indices. The fact that we are dealing with su{2) rather 
than su(l, 1) implies that the two-metric q a b : = e l a e\ has Euclidean signature. Moreover, 
let A\ be an su(2) connection and define the field Ef := e ab e l b where e a b is the metric- 
independent totally skew tensor of density weight —1. Then it turns out that the pair 
{A l al E?) is a canonical one for the Hamiltonian formulation of 2+1 gravity based on the 

Einstein Hilbert action S = Sm^ x \J\ det(g)\R^ where g is the three-metric and 

its scalar curvature. In other words, Ef is the momentum conjugate to A l a so that the 

symplectic structure is given by 

{K{x),E b {y)} = 5 b a 5)5{x,y). (2.1) 

The Hamiltonian of the theory is a linear combination of constraints, / d 2 x(A l Gi + N l Ci) 
for some Lagrange multipiers h % ,N % where 

d := D a E°; = d a E°; + eij k A{El : Gauss constraint, 

d := ^e ab F l ab : Curvature constraint (2.2) 

where F a & denotes the curvature of A a . The Gauss constraint appears also in 3+1 gravity, 
however, the curvature constraint is completely different from the constraints that govern 
3+1 gravity, ||. The equivalent of Ci in 3+1 gravity are two types of constraints, one of 
them, V a , generates diffeomorphisms, the other one, H, generates dynamics. The curva- 
ture constraints Ci on the other hand do not generate any such gauge transformations, 
in fact, the connection Poisson-commutes with Cj and shows that it is a Dirac observable 
with respect to Cj. The constraint Ci = imposes that the connection should be flat and 
thus the classically reduced phase space becomes the cotangent bundle over the moduli 
space of flat su(2) connections which is finite-dimensional. 

It is obvious that the quantization of the model as defined by fl2.2|) will not give too much 
insight into the 3+1 situation. In the following we will reformulate ( |2.2| ) in such a way 
that it brings us in connection with 3+1 gravity. 

It will turn out that the following compound field, called the degeneracy vector, for reasons 
that will become obvious soon 

E l := X -e^ k e ah E a 3 E\ (2.3) 
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is a crucial one. Let us compute the square of this density of weight one 



= lq bd [E b 3 Ef] = l -q bd qa C e ba e dc 

= det(g), (2.4) 

that is, the two-metric is degenerate if and only if E± = is identically zero. We also see 
that det(g) is manifestly non- negative. Notice that E^E % = 0. 

Whenever the degeneracy vector is non-vanishing we can perform the following non- 
singular transformation (N l ) <-> (N a ,N) for a vector field N a , called the shift, and a 
scalar function N, called the lapse : 

jTji E^E a E^ 

N l = N a e ab E b + N ^N a = e tjk —/-N k , N = . (2.5) 

/det(g) det(g) Jdet(g) 



Notice that formula fl2.5|) respects that N l ,N a ,N have density weight zero. Using ( |2.5|) , 
we can now write the curvature constraint in the form 

N l d = N a V a + NH where 

V a := F^ b E b : Diffeormorphism constraint 

1 E a E^ 

H := -Fg b ejjk 3 '■ Hamiltonian constraint . (2-6) 

2 ^det(g) 

Apart from the fact that we are in two rather than three space dimensions these are pre- 
cisely the constraints of Euclidean 3+1 gravity 0. Since the 3+1 Euclidean Hamiltonian 
constraint operator plays a key role in the quantization of the 3+1 Lorentzian Hamilto- 
nian constraint we claim that the set of constraints ( |2.6| ) bring us in maximal contact 
with the 3+1 theory. 



Notice that unlike in Q], we have a factor of l/ydet(g) in the definition of the Hamiltonian 
constraint. This difference has two important consequences : 

1) Classical : 

The denominator in H = F i E i / v /det(g) where F { := e ab F^ b /2 (or F* b = e^F*) 

blows up as E % vanishes. Since the limit lim^ E/\ \E\ \ depends on the details of 
the limiting procedure we must exclude degenerate metrics classically. This is in 
contrast to || where the authors exploit the possible classical degeneracy of the 
metric when one discards the denominator to demonstrate that one has already an 
infinite number of degrees of freedom at the classical level (notice, however, that 
their solutions, where F % or E l become null, do not apply since we are dealing with 
su(2)). 

2) Quantum : 



It is by now known that one of the reasons for why H := Jdet(q)H suffers from 



1 In fact, once one has densely defined the Euclidean operator in 3+1 dimensions the fact that the 
Lorentzian operator in 3+1 is densely defined is a simple Corollary l|. 
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huge problems upon quantizing it is due to the fact that H has density weight two 
rather than one. As argued in || , only densities of weight one have a chance to be 
promoted into densely defined, covariant operators. This is why we must keep the 



denominator l/y'det(g) in ( |2.6| ) at the quantum level. 

Just like in 3+1 gravity we wish to work in a connection representation, that is, states are 
going to be functions of connections. Then an immediate problem with H is that one has 



to give a meaning to the denominator l/ydet(g). In || that was achieved for Lorentzian 
3+1 gravity by noting that the denominator could be absorbed into a Poisson bracket 
with respect to a functional V of q a b. The idea was then to use the quantization rule that 
Poisson brackets should be replaced by commutators times 1 / (iti) and to replace V by an 
appropriate operator V. Such an operator indeed exists and it is densely defined. 

Is a similar trick also available for 2+1 gravity ? At first sight the answer seems to be 
in the negative because the underlying reason for why such a trick worked for 3+1 gravity 
was that the co-triad e l a , the precise analogue of the degeneracy vector E\ considered as a 
function of Ef was integrable, the generating functional being given by the total volume 
V of S. In other words, we had 



1 Tpb Tpc AT/ r 

t^k -^k = |^ with V : = / d 3 xjfet(q) (2.7) 

2 Jdet(q) 5E i 7s 



However, if we take over the definition of V (with d 3 x replaced by d 2 x) then we find 
instead 

SV q ah E\ 
/det(g) 

Thus, there seems not such a trick available in the 2+1 case. However, it is a matter of 
straightforward computation to verify that indeed 

E* = ±e a % k {Ai,V}{A k b ,V} (2.9) 



which does not seem to help much because what we need is E l /^det(q) rather than E % 
itself. 

The new input needed here as compared to the 3+1 case is as follows : Notice that if we 



could replace ydet(g) by V then we could absorb it into the Poisson brackets by using 
the identity 

{Ai,V}{A k b ,V} 



y 4{Ai,VV}{A«,VV} 

As we will see, V can be promoted, just as in the 3+1 case, into a densely defined positive 
semi- definite operator. Therefore its square root exists and it would follow that the last 
equation with Poisson brackets replaced by commutators would make sense as an operator. 
In the next section we will define a Hilbert space and the corresponding operator. 



What remains is to justify the replacement of ydet(g) by V. That this is possible we will 
show in the section after the following. It happens because the Poisson bracket gives a 
local quantity and therefore we may actually replace V by V(x, e) in 

{Ai(x), V} = {Ai(x),V(x,e)} where V(B) = [ d 2 x^fdet(q) 
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is the volume of a compact region B and V(x, e) is the volume of an arbitrarily small open 
neighbourhood of the point x, the smallness governed by e. It is then easy to see that 
lim e ^o V(x, e)/e 2 = Jdet(q) (x). Now, in the quantum theory we are going to point split 
the quantity H and we will use a regularized 5 distribution with point split parameter e. 
As we will see, that parameter can be absorbed into V(x, e) to serve as a replacement for 
\J (det(q)(x). The details are displayed in the following sections. 

3 Quantum Theory and Volume operator 

In this section we will review the definition of a Hilbert space for diffeomorphism invariant 
theories of connections || . This will be our kinematical framework. On that Hilbert space 
we are going to construct a 2+1 volume operator which turns out to be actually more 
complicated than the one for the 3+1 theory 0. [|. 

3.1 Quantum kinematics 

In what follows we give an extract from [0, [F]. The reader interested in the details is 
urged to study those papers. 

We will denote by 7 a finite piecewise analytic graph in E. That is, we have analytic 
edges e which are joined in vertices v. We subdivide each edge into two parts and equip 
each part with an orientation that is outgoing from the vertex (the point where these two 
parts meet is a point of analyticity and therefore not a vertex of a graph, thus each edge 
from now on can be viewed to be incident at precisely one vertex). Given an su(2) con- 
nection A l a on E we can compute its holonomy (or path-ordered exponential) h e (A) along 
an edge e of the graph. Recall that all representations of SU(2) are completely reducible 
and that the (equivalence class of equivalent) irreducible ones can be characterized by a 
half integral non- negative number j, the spin of the representation. We will denote the 
matrix elements of the j-representation at g 6 SU(2) by ftj(g)- 

Consider now a vertex v of the graph and the edges e%, .., e n incident at v, that is, the 
graph has valence n. Under a gauge transformation g at v the holonomy transforms as 
h ei — > gh ei , i — 1, ..,n. Now consider the transformation of the following function 

®i=^k{K) -> [®u*M] ■ [®r=i^ej] • 

We are interested in making this function gauge invariant at v. To that end we or- 
thogonally decompose the tensor product of the Hj^g) into irreducibles and look for the 
independent singlets in that decomposition. There is an orthogonal projector c v on each 
of these singlets, we say that it is compatible with the spins ji, ..,j n , and so we can make 
our function gauge invariant at v by contracting : c v ■ ®i=i'n'j i {h ei ). 

If we do that for each vertex we obtain a completely gauge invariant function called a 
spin-network function. Thus a spin-network function is labelled by a graph 7, a colouring 
of its edges e with a spin j e and a dressing of each vertex v with a gauge- invariant projector 
c v . If we denote by E(j), V(j) the set of edges and vertices of 7 respectively then we use 
the shorthand notation 

T^g where j := {j e } eejB(7 ), c := {c v } veVM , 

for that spin-network function. 
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The Hilbert space Ji that we are going to use for gauge invariant functions of con- 
nections is most easily described by saying that the set of all spin-network functions is 
a complete orthonormal basis of 7i (so each spin-network function comes with a specific 
finite normalization factor). Notice that therefore TC is not separable. Another charac- 
terization of 7i which is very useful is to display it as a certain L2 space. To that end, 
consider the finite linear combinations $ of spin-network functions. $ can be turned 
into an Abelian C* algebra by saying that involution is just complex conjugation and by 
completing it with respect to the sup— norm over the space A/Q of smooth connections 
modulo gauge transformations. That C* algebra is isometric isomorphic by standard 
Gel'fand techniques to the C* algebra of continuous functions C(A/Q) where A/Q is the 
set of all homomorphisms from the original algebra into the complex numbers. The space 
A/Q, as the notation suggests, is a certain extension of A/Q and will be called the set 
of distributional connections. Indeed, it is the maximal extension such that (the Gel'fand 
transform of the) spin-network functions are continuous. By standard results, the result- 
ing topology is such that A/Q is a compact Hausdorff space and as such positive linear 
functionals T on C(A/Q) are in one to one correspondence with regular Borel measures 
fi on A/Q via r(/) =: //(/) = fjjgd/if. 

Now the measure fiQ underlying 7i is completely characterized by the integral of spin- 
network functions and is given by HoiT^j^ = 1 if T 7 j s = 1 and otherwise. So we have 
Tt = L 2 (A/Q,dfi ) and spin-network functions play the same role for /i that Hermite 
functions play for Gaussian measures. 

In the sequel we will topologize the space $ of finite linear combinations of spin- 
network functions in a different way and we will call $ henceforth the space of cylindrical 
functions. A function / 7 is said to be cylindrical with respect to a graph 7 if it can be 
written as a finite linear combination of spin-network functions on that 7. The norm of 
/ 7 will be the L\ norm ||/ 7 ||i = J2i \ < Ti, f > | which equips $ with the structure of 
a topological vector space. The distributional dual $' is the set of all continuous linear 
functionals on $. Certainly every element of Ti is an element of $' by the Schwarz 
inequality and every element of <3> is trivially an element of Ti. Thus we have the inclusion 
$ C W C $' (this is not a Gel'fand triple in the strict sense because the topology on $ is 
not nuclear). 

This furnishes the quantum kinematics. Notice that we can take over the results 
from || without change concerning the Diffeomorphism constraint : given an analyticity 
preserving diffeomorphism ip we have a unitary operator on 7i which acts on a function 
cylindrical with respect to a graph as {7(<^)/ 7 = f v M, that is, the diffeomorphism group 
Diff(S) is unitarily represented. This implies that one can group average with respect to 
the diffeomorphism group as in ||. We will return to this point in section 6. 

3.2 The 2+1 volume operator 

The plan of this subsection is as follows : 

Since E l (x) is a density of weight one, it makes sense that it will give rise to a well-defined 
and diffeomorphism-covariantly defined operator valued distribution. In a second step we 
will point-split det(g) = E l E l and take the square root of the resulting operator. Again, 
since y / det(g) is a density of weight one, it can be turned into a well-defined operator- 
valued distribution even in regulated form and the limit as the regulator is removed exists. 
Let us then begin with E l . Let as in the previous section f 1 denote a function cylindrical 
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with respect to a graph 7 and denote by E(j) its set of edges. Edges are, by suitably sub- 
dividing them into two halves, in the sequel always supposed to be oriented as outgoing 
at a vertex. We will compute the action of various operators first on functions of smooth 
connections and then extend the end result to all of A/Q. 

Let 5g(x,y) = 5 tl (x 1 ,y 1 )5 t2 (x 2 ,y 2 ) be any two-parameter family of smooth functions of 
compact support such that lim ei)£2 _ > o is d 2 yS e (x, y)f(y) = f(x) for any, say smooth, func- 
tion on E where e = (ei, e 2 ) parametrizes the size of the support. Consider the point-split 
operator 

Ei?{x) := \e ab e^ k £ d 2 y £ d 2 z5,(x, y)5?(x, z)E«(y)E b k (z) (3.1) 
and apply it to / 7 . Notice that upon replacing Ef{x) = —ih5/5A l a (x) 

E?(x)f, = -th Yl f'dtS&eitye^XKt)^ (3.2) 

eg£( 7 ) J0 

where X l e {t) = tv([h e (0, t)Tih e (t, l)] T d/dh e (0, 1)), h e (a, b) is the holonomy from parameter 
value a to b and r, are generators of su(2) with structure constants e^. We also need the 
quantity 

X l J\s,t) = tT([h e (0 1 s)nh e (s,t)T j h e (t 1 l)] T d/dh l ,(0 1 i)) for s < t (modulo 1). Then it is 
easy to see that 

= -?i 2 ^e ab e ijk [ d 2 y f d 2 z5 € -{x,y)5?{x,z) x 

dt5(y,e(t))e a (t) dt'5(z,e'(t'))e' b (t') x 
Jo 

x \X k e ,{t')Xi{t) + 8 e , e ,{8{t,t')Xi k (t,t') + 6{t',t)Xi k {t',t)}}f^ 

= -h 2 -t ah t ijk V f 1 dte a (t) f 1 dt'e' b (t')5g(x,e(t))8?(x,e'{t')) x 

2 e,e'ei?(V J ° 

x [X k (t')Xl(t) + 5 e , e ,{0(t,t')Xl k (t,t') + d(t',t)Xi k (t',t)}]f, (3.3) 
where 9(s,t) — 1 if s < t and otherwise. 

We are now interested in the limit e — > and proceed similar as in |J. We must adapt 
the regularization to each pair e, e' to get a well-defined result. 

1) Case e = e'. 

If x does not lie on e then for sufficiently small e we must get <5?(x, e(t)) = for 
any t G [0, 1]. Thus in the limit we get a non- vanishing contribution if and only 
if there exists a value t x e [0, 1] such that e(t x ) = x (there is at most one such 
value t x because edges are not self-intersecting). Since e is nowhere vanishing we 
must have e 1 ^) 7^ (switch 1 <-> 2 if necessary). We send e%, e[ — > and find that 
6g(x,e(t)) —> 5 t2 (x 2 ) e 2 (t))5(t — t x ) /\e l (t x )\ and similar for 5p(x, e'(t')). Inserting this 
into ( |3.3| ) we find that there is no contribution for e = e' because of the two zeroes 
= €abe a (t x )e b (t x ) and = eijk[X l J{t x , t x ) +X^{t x , t x )]. Notice that it was crucial to 
have e 2 , e' 2 still finite as otherwise the appearing S e2 (0)8 e > (0) would be meaningless. 
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2) Case e ^ e'. 

If again x does not lie on both e, e' then by choosing e, e* sufficiently small we 
must get zero. Therefore e, e' must intersect and as we have divided edges into two 
halves they can intersect at most in their common starting point corresponding to 
t = t' = which is thus a vertex v of the graph 7. 

A) Subcase 

Consider first the case that e, e' have co-linear tangents at t — and let us 
assume that e 1 (0),e' 1 (0) 7^ (switch 1 <-> 2 if necessary). Then we first send 
£1, e' x — > which results in 

2 „2/-xxr 2 „/2/./^ <*(*)<*(*') 



^,e(t))^(x,e'(t')) - M^(t))5 4 (xVl*')) 



|e 1 (0)e /1 (0)| 



and thus performing the two t integrals we get zero as above because = 
£afee a (0)e /fe (0) by assumption. 

B) Subcase 

We are left with the case that the tangents of e, e' are linearly independent at 
x = v . We replace 8?(x, e(t))5?(x, e'(t')) by 6?(e'(t'), e(t))5?(x, v) and send first 
e — > 0. Then 

He'd') e(t)) - mm 

and we can perform the integral. Since we are integrating over a square [0, l] 2 
and the two-dimensional delta-distribution is supported at a corner we pick up 
a factor of 1/4 upon setting t = t' = and dropping the integral. At last we 
send e f — > 0. 

Summarizing, we find (^(7) denotes the set of vertices of 7) 

E\x)f, = -^- £ 5(x,v) £ sgn(e,e')e^X e fc ,/ 7 (3.4) 

u€V(7) e,e'G-E(7),ene'=i; 

where X* := ^'(0) is easily recognized as the right invariant vector field on SU(2) eval- 
uated at g = h e (0, 1) and sgn(e, e') is the sign of e a be a (0)e' b (0) and so is an orientation 
factor. This furnishes the definition of the operator corresponding to the degeneracy 
vector. 

We now will define the volume operator for any compact region BcE. Our first task 
is to define an operator corresponding to det(g) and then to take its square root. Since 
det(g) is a density of weight two we expect this to be quite singular, in fact the naive 
definition det(q)(x) := E t (x)E l (x) does not make any sense given the expression (|3.4|) 
which involves a factor of 5(x,v). Thus we are lead to point-split the two degeneracy 
vector operators and to hope that 1) the regulated operator is positive so that it makes 
sense to take its square root and 2) that one can remove the regulator from the square 
root. Let us then define similar as above 

detfg) e> -,(x) := fj 2 y J d 2 z5 e -(x,y)5 e ,(x, z)E\y)E\z) (3.5) 
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and apply it to a function cylindrical with respect to a graph 7. Given ( |3.4j ) the result is 
easily seen to be 

det(g)g. ? (a;)/ 7 
= g4 E 8g(x,v)5?(x,v') x 

t),»'6V(7) 

x £ £ sgn(e 1 ,e 2 )sgn(ei,4) x 

ei,e 2 e-B(7),eine2=f e'^e'^E^fi^He'^v' 

x [e^XH eimn ^^lA- ( 3 - 6 ) 

We now will accomplish both hopes 1), 2) stated above by appropriately choosing the 
regulators. 

1) Choose e =: e*, then we are able to display ( |3.5|) as a square of an operator 

det(g)- ? (x)/ 7 

= {y E E sgnle^e^^X^]} 2 /,. (3.7) 

Since X % X° e commute for e\ 7^ e2 and because %X\ is essentially self-adjoint with range in 
its domain, so is X\ Jf| 2 and therefore the whole operator corresponding to one factor in 
(p.7|). Thus, ( |3.7| ) is a square of essentially self-adjoint operators with range in its domain 
and so it is positive semi-definite. Therefore its square root is well defined. 

2) Choose e small enough such that 5g(x, v)5g(x, v') = 5 v y[5g(x,v)} 2 , that is, given 7,2; 
we must choose e so small that for v ^ v' not both of them can be in the support of the 
function 5?(x, .) which is always possible. Then we may write (|3.7|) as 



det(g) e> -(x)/ 7 

^ £ fefo*)] a { E sgn( ei ,e 2 )^^X e fc 2 ]} 2 / 7 , (3.f 



USV(7) ei,e2G-E(7),eine2=f 



take its square root and define this to be the regulated operator corresponding to Jdet(q) 



y / det(^) e ,(x)/ 7 := ^t^y^f, • (3.9) 

In considering the limit e — > notice that for small enough e at most one vertex of 7 
lies in the support of 8g(x, .). Therefore we can take the sum over vertices and the factor 
[Sg(x, v)} 2 out of the square root and find that 



2 



det( g y*)/ 7 = y £ 8&,v\K E sgn(e 1 ,e 2 )[6^X| 1 X4]}V 7 . (3.10) 

v€V(y) y ei,e2G-B(7),eine2=u 

But now the limit e — > is trivial to take, we finally find that 

— h 2 



'det(g)(x)/ 7 = - £ 6(x,v) { E sgn( ei , e 2 )[e^Xg 1 X e fc J} 2 i 7 (3.11) 

«6V(7) V ei,e 2 e-E(7),eine2=D 
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or in integrated form 




V(B)f,:=ll d 2 x^Jdet(q)(x)]l 

B 



{ £ sgn(e 1 ,e 2 )[e^X| 1 X^]}V 7 . (3.12) 

ei,e2£-E(7),eine2=f 



Formula (|3.12j ) motivates to introduce the "volume operator at a point" V v : For each 
integer n > 2 define {[i>, Ji]} to be the set of germs of n analytical edges incident at v 
(a germ of an analytical edge at a point v is a complete set of analytical data available 
at v that are necessary to reconstruct it, that is, essentially the coefficients of its Taylor 
series). For a germ e n := (ei, ..,e n ) G {[t>,n]} define 



V e J"i e ./Ge 



where the right invariant vector field X l e (g) = X\(gKfih € SU(2), due to right invariance, 
depends really only on the germ of the edge e because it acts on a function in the same 
way no matter how "short" the segment of e is on which that function actually depends, 
as long as that segment starts at v = e(0). In particular all XI, e incident at v, commute 
as long as their germs are different. Then 

oo 

V(B) = Y,V V where V v = £ £ V Sn . (3.13) 

v&B n=2 e n e{[v,n]} 

We see that V(B) is a densely defined, essentially self-adjoint, positive semi-definite op- 
erator on 7i for each bounded region BcE. Its most interesting property is that it acts 
non-trivially only at vertices of the graph underlying a cylindrical function, moreover, 
that vertex has to be such that at least two edges incident at it have linearly indepen- 
dent tangents there. This is in complete analogy with the volume operator of the 3+1 
theory just that we need to replace everywhere valence three by valence two. Unlike 
the the three-dimensional volume operator, however, its two-dimensional "brother" does 
not vanish at two-valent and three-valent vertices at all as long as there are at least two 
edges with linearly independent tangents at the vertex under consideration. As we will 
see in the appendix, the two-dimensional volume operator is even positive definite on 
gauge invariant functions with two-and three valent vertices while the three-dimensional 
volume operator annihilates such functions identically. This is to be expected because 
by inspection of ( |3.13j ) the principal symbol of that operator is non-singular on two-and 
three valent vertices while in the three-dimensional case it is singular. 
The fact that the volume operator acts only at vertices of the graph will enable us to 
take the infra-red limit in case we are dealing with asymptotically flat topologies and also 
ensure that the ultra-violet limit exists. Thus, the volume operator acts both as an IR 



and as an UV dynamical regulator, a point of view emphasized in flQfl . 



Remark : 

Notice that q a b = {A l a , V}{A* b , V} just as in the three-dimensional case. This observation 



lead in the three-dimensional case to the construction of a length operator [11]. The 



only crucial property that was necessary to construct this operator was that the volume 
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operator acts only at vertices. Since that is true for the two-dimensional operator as well 
we can therefore take over all the results and formulae from [11 J to the two-dimensional 
case, except for the obvious differences which are due to different dimension and algebraic 
expressions in terms of right invariant vector fields of the volume operators. In particular, 
although the eigenvalues of the length operators are certainly different, qualitatively the 
spectrum is still discrete, the operator is positive semi-definite and essentially self-adjoint 
and the length of a curve as measured by a spin-network state is different from zero only 
if at least one edge of the graph crosses the curve, though not necessarily in a vertex. 
Thus we automatically have a two-dimensional length operator as well. 
The fact that the two-dimensional length operator is less degenerate than the three- 
dimensional one can be traced back to the observation that what is length in two dimen- 
sions is what is area in three dimensions. 



4 Regular izat ion 

This section is divided into three parts : in the first part we will derive a regulated 
Wheeler-DeWitt operator. The regularization consists in a triangulation of £ which is 
kept arbitrary at this stage. In the next part we will specify the properties that we wish 
to impose on the triangulation and then make a particular choice which satisfies those 
properties. Finally in the last part we complete the regularization by employing that 
triangulation and take the continuum limit which then equips us with a densely defined 
family of operators, one for each graph. 

The presentation will be kept largely parallel to the one in H] in order to fasciliate com- 
parison. 



4.1 Derivation of the regulated operator 

We wish to define an operator corresponding to 

H(N) := [ d 2 xNF i ^= 
■fc v /det(^) 

- -/^ M ) (4.1) 

where we have used that tr(TjTjTfc) = — 6^/2 and (|2.9|) . Following the idea outlined in 
section 2 consider now a point splitting of the above expression as follows : Let e be a 
small number and Xe{ x , u) '■= $(§ — l^ 1 ~ l/ 1 |)^(f — \ x<2 ~ U 2 \) where 6{t) — 1 if t > and 
otherwise, that is, Xt is the characteristic function of square of coordinate volume e 2 . 
Moreover, it is just true that {A l a (x), V} = {A l a (x), V(x, e)} where 

V{x,e) := / d 2 yxe(x,y)\/det(q)(y) 
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is the volume of the square around x as measured by qw,. Notice that trivially lim e ^ V(x, e ) / g2 
det(q)(x). Therefore we have the identity (we write the density F % as a 2-form) 

{A(y),V}A{A(y),V}, 



limy N(x)tr(F(x) J Xe(x,y) 
lim J N(x)tr(F(x) J 



Xe(x,y) {A(y),V}A{A(y),V} 
e 2 V(y,e)/e* } 

N(x)tr(F(x) [ lhm X^ ] iMlXlAiMlXl ) 
= ~H(N), (4.2) 

that is, the point splitting singularity 1/e 2 was absorbed into V(y,e). The limit identity 
( |4.2| ) motivates to define a point split expression 

H,(N) := -/*<«)*<*■<,) 



4 / N(x)tr(F(x) f x ,(x, y){A(y), Jv(y,c)} A ^Vd/.e)}), (4.3) 



that is, the simple formula {., yV(y, e)} = {., V(y, e)}/(2V(y, e)) enabled us to bring the 
volume functional from the denominator into the nominator, of course, inside the Poisson 
bracket. 

The idea is now to replace Poisson brackets by commutators and the volume functional 
by the volume operator and then take the limit e — > 0. In order to do that we must first 
write (O) in such quantities on which the volume operator knows how to act. Since, as 



obvious from the previous section, it only knows how to act on functions of holonomies 
along edges we must replace the connection field A l a in ( [4.3[ ) by holonomies. We are thus 
forced to introduce a triangulation of S. 

Denote by A a solid triangle. Single out one of the corners of the triangle and call it 
v (A), the basepoint of A. At v(A) there are incident two edges Si(A), s 2 (A) of dA which 
we equip with outgoing orientation, that is, they start at v(A). We fix the labelling as 
follows : let s be the analytic extension of si(A) and si(A) the half of s starting at v(A) 
but not including sj(A) — {v} with outgoing orientation at v(A). Let U be a sufficiently 
small neighbourhood of v(A) which is split into two halves by s. Define the upper half 
U + of U to be that half of U which one intersects as one turns si(A) counterclockwise 
into si(A). Now we require that there exists U such that U fl 52(A) = U + H 52(A), that 
is, 82(A) intersects the upper half of U. 

Definition 4.1 Two analytical edges ei,e 2 incident and outgoing at v = e\ fl e2 will be 
said to be right oriented iff there exists a neighbourhood U of v, its upper half U + being 
defined by e\, such that e 2 intersects U + . 

This prescription is obviously diffeomorphism invariant. Notice that we did not, as it is 
usually done for triangulations, require that the tangents of the edges bounding A must 
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have linearly independent tangents at their intersection. If they are linearly independent 
then our prescription is equivalent to saying that e a i>s± (0)53(0) > 

Finally, let a(A) denote the remaining edge of <9(A), called the arc of A, whose orientation 
we fix by requiring that it runs from the endpoint of si(A) to the endpoint of 52(A). Then 
dA = 012(A) = si(A) o a (A) o s 2 (A)~ 1 is called the loop of A based at v(A). We define 
also a 2 i(A) := ai 2 (A) _1 . 

Let us now write the integral over E x E in ( f4.3| ) as a double sum of integrals over 
A x A' where A, A' are triangles of some triangulation T of S 




(4.4) 

The purpose of the notation just introduced is that we may approximate, for sufficiently 
fine triangulation, each of the integrals by a function of holonomies as follows : Let 5 be 
a small parameter and s«(A) be the image of [0,5] under the path Sj(A, t). Then, using 
smoothness of the connection we find 

N(v(A r )) X e(v(A'),y)e^h aij{A/) 

= 25 2 N(v(A'))xM*), V) * ?(A '' °y (A/? 0) F ab {v{A')) + 0(6") 
= 2/ Xe(x,y)N(x)F(x) + o(5 3 ) and 

J A' 

X,(x,v(A))e^h St{A) {h^ A) , ^V(v(A),e)}h S]{A) {h^ {A) , y/v(v(A),e)} 
= X e(x, viA^s^A, 0)4(A, 0){A a (v(A)), \Jv(v(A), e)}{A b (v(A)), \Jv(v(A), e)}) 
= 2 J X e{x, y){A(y), ^JV(y,e)} A {A(y), ^Jv(y,e)} + o(5 3 ) (4.5) 

since the area of A is approximately 5 2 e a bs1(A, 0)s|(A, 0)/2 so that both integrals are of 
order 5 2 provided that the tangents of <9(A) at v(A) are linearly independent. Thus, up 
to an error of order 5 2 which vanishes in the limit as the we remove the triangulation we 
may substitute ( |4.4|) by 

H T , e (N) = -2 e^e H iV(t;(A'))Xe(^(A / ),t;(A))x 

A,A'eT 

x tr(/i ai . (A0 /i Sfc (^{^(A)' V /K ( t; ( A )' e )}^(A){^(A) ) \fv(v(A),e)}) . 

(4.6) 

The result ( |4.6| ) is still purely classical and becomes H(N) when taking 

1) first the continuum limit (that is, refining the triangulation ad infinitum) and 

2) taking e — > on smooth connections A\ and smooth momenta 

A second way to guide the limit and that leads to H(N) is by "synchronizing" e ~ 5 and 
to take 5 — > as follows : for each A define 

e(A) := v /|e afe sf(A,0)^(A,0)|5, 

replace for each A' : 

l) X e(v(A),v(A')) by Xe(A')(v(A),v(A')) and 
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2) V(v(A'),e)byV(v(A'),e(A')) 

and then take 5 — > 0. Notice that this corresponds to introducing e(y) = p(y)S instead of 

e in ( |4.2| ) where p(y) is an almost nowhere (with respect to d 2 x) vanishing function such 

that p(v(A))5 = e(A). Clearly p must be almost nowhere vanishing as otherwise we do 

not get a 5 distribution in the limit 5 — > 0. Notice that the set of v (A')'s has d 2 x measure 

zero so that a vanishing p(v(A') is not worrysome. 

It will be this latter limit which is meaningful in the quantum theory. 

We have managed to write H(N) in terms of holonomies up to an error which vanishes 

in either of the limits that we have indicated. 

The next step is to turn ( }4.6|) into a quantum operator. This now just consists in 
replacing ^(^(A^e) by V(v(A),e) and Poisson brackets by commutators times l/{ih) 
because we work in a connection representation. The result is 

H T>e (N) = ^ E ^e kl N(v(A f ))xM^'),v(A))x 

x tr(h ai . {A/) h Sk{A) [h; k \ A) , <Jv(v{A), e)]^ ( A)[^; A) , <Jv(v(A), e)]) . (4.7) 

We wish to show that ( |4.7| ) is densely defined in the limit e — ► no matter how we 
choose the triangulation T, as long as it is finite, thereby showing that the regulator 
e can be removed without encountering any singularity. Thus, we prescribe the e — > 
limit before taking the limit of infinitely fine triangulation (continuum limit) and therefore 
have interchanged the order of limits as compared to the classical theory. However, as we 
will show shortly, one arrives at the same result when synchronizing e ~ 5 and taking 8 
sufficiently small but finite for the moment being which corresponds to the second way to 
guide the classical limit indicated above and therefore interchanging the limits is allowed. 
For that purpose let / 7 be a function which is cylindrical with respect to a graph. Consider 
first some triangle A which does not intersect 7 at all. Then it is easy to see that 



KUJvmau)])/, 







The reason for this is that the graphs 7 and 7 U si(A) then do not have any two-valent 
vertex in the box around v(A) parametrized by e other than the vertices of 7 themselves. 
Thus the volume operator does not act on hj^^ and the commutator vanishes. It follows 
that only tetrahedra which intersect the graph contribute in ( f4.7|) . So let 7 H A 7^ 0. For 
the same reason as above we find a non-zero contribution only if Si(A) or 52(A) intersect 
7, that 012(A) alone intersects 7 is not sufficient. Moreover, still for the same reason, 
if Si (A) intersects 7 but not in the starting point of Sj(A) then we still get zero upon 
choosing e sufficiently small so that the intersection point p lies outside the support of 
the characteristic function, that is, Xe(v(A),p) = 0. Thus a triangle A contributes to 
( |4.7| ) if and only if v (A) e 7. But if that is true then we may replace V(v(A), e) by the 
operator K(A) defined in ( |3.13|) and so the e-dependence of V(v(A),e) has dropped out. 



The remaining e-dependence now just rests in the function Xe(v (A'), v(A)). Now, since 
we let e — > first, at finite triangulation, we conclude altogether that the unrestricted 
double sum over triangles in ( [4.7] ) collapses to a double sum over triangles subject to the 
condition that their basepoints coincide and lies on the graph. In formulae 

i^ T (iV)/ 7 :=lim^ Tie (iV)/ 7 
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A,A'67>:=u(A)=u(A')e7 
X ^(^(AO^^A)^;^), V^]^(A)[^(A)> "V^D/y ( 4 - 8 ) 

which displays Ht(N) as a densely defined operator which does not suffer from any sin- 
gularities because at finite triangulation there are only a finite number of terms involved 
in ( |4.8|) , even if £ is not compact. 

Notice that in the e — > limit we have recovered a gauge invariant operator as we should. 



Let us now show that one arrives at the same result by synchronizing e « 5 as above 
and taking 5 sufficiently small but still finite : Namely, by choosing e(A') as above we 
have arranged that only the starting points of the Sj(A') are covered by the e(A')-box 
around v(A') that underlies the definition of V(v(A'), e(5')). This implies first of all that 
we need to sum only over v(A) = v(A'). Next, as we will be forced to adapt the trian- 
gulation to the graph anyway, we can arrange that the A intersect 7 only either in whole 

edges or in vertices of A. If that is the case, then it follows that [h Si (A), \fv( v (A), e(A))]/ 7 
is non-vanishing only if s$(A) intersects 7 in v(A) because the end-point is not covered 
by the e(A)-box and if s$(A) is contained in 7 but does not start in a vertex of 7 then 
the commutator vanishes due to the properties of the volume operator. This is enough to 
see that we arrive at fl4.8| ) again. 

In either way of taking the limit we are now left with taking the continuum limit 5 — > of 
refining the triangulation ad infinitum which we denote as T — > 00. Certainly that limit 
depends largely on the choice of the limit T — > 00. For instance, if we are not careful and 
refine T in such a way that the number of basepoints of triangles that intersect 7 diverges 
we will not get a densely defined operator. We see that we must choose T according to 7 
so that we get actually a family of operators 

H JjT (N) = H T{J) (N) 

where T(j) is a triangulation adapted to 7 together with a well-defined refinement proce- 
dure T — > 00. We will propose such a T(j) in the next subsection guided by some physical 
principles. It will then be our task to verify that the family (H 1j t(N)) still defines a linear 
operator. 



4.2 Choice of the triangulation 

So far everything what we said was in complete analogy with the three-dimensional case 
except that there we did not even need a point-splitting. In particular, (|4.8|) is the 
precise counterpart of the three-dimensional Euclidean Hamiltonian constraint operator. 

What is different now is that in the 3+1 case the volume operator was much more 
degenerate than in the 2+1 case, a result of which was that a basepoint of a simplex had 
to coincide with a vertex of the graph in order to contribute without further specification 
of the triangulation. Therefore, it was sufficient to adapt the triangulation to the graph in 
such a way that, among other things, the number of simplices intersecting a vertex stays 
constant as one refines the triangulation in order to arrive at well-defined continuum 
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limit. In the 2+1 case that is not true any longer and one must worry about the number 
of triangles intersecting the graph 7 off the vertices of 7. 

Let us adopt the physical principles listed in |J which should guide one of how to 
choose the triangulation. In brief, they were : 

1) The amount of ambiguity arising from the choice of the triangulation should be kept 
to a minimum. 

2) The resulting operator should be non-trivial and not annihilate every state. 

3) The choice of the contributing A should be diffeomorphism covariant as to interact 
well with the diffeomorphism invariance of the theory. 

4) The choice of the A should be canonical and not single out one part of the graph 
as compared to the other or one graph as compared to another. 

5) The family of operators H^(N) should define a linear operator H(N) (cylindrical 
consistency) . 

6) The resulting operator H y<T (N) should be densely defined with a well-defined con- 
tinuum limit. That is, if $ G $' is a diffeomorphism invariant distribution and / 7 a 
function cylindrical with respect to a graph 7 then 

lim V(H^ T (N)f 7 ) =: 9{H(N)fJ) 

1 — >oo 

exists. The fact that ^ is diffeomorphism invariant is because we actually want to 
define H(N) on solutions to the diffeomorphism constraint which turn out to be 
distributions and so the above limit is the precise sense in which H(N) is defined 
on distributions. 

7) The operator H(N) should be free of anomalies, that is, 

*([H{M),H(N))<f>) = 
for each <p G $ and every diffeomorphism invariant G 

Since we wish to obtain a densely defined operator no matter how fine the triangulation 
while keeping the extra structure coming from the triangulation to a minimum we are 
naturally lead to impose that the triangles that intersect 7 in its basepoint must be 
constant in number. There are only two diffeomorphism invariantly different possibilities 
: either v(A) is a vertex of 7 or it lies on an edge of 7 between its endpoints. Since we 
want to get a non- vanishing operator one of the two or both scenarios should happen. 

Suppose first then that v(A) is an interiour point of an edge e. Then there is no 
natural way how to choose the triangle A itself : the only structure available is the edge 
e and one may therefore choose one of Sj(A), say si(A), to lie entirely in e. But then 
52(A) should certainly not lie in e otherwise v(A) would be a vertex of 7 U dA with only 
co-linear tangents of edges incident at it and the volume operator V v (a) would vanish. 
Thus there is at least a huge ambiguity in how to choose 52(A). 

If, on the other hand, v(A) is a vertex of the graph then there are at least two edges 
e, e' of 7 incident at it and now it is a natural choice to assume that Sj(A) coincide with 
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segments of e,e'. 

In conclusion, guided by the principle of introducing as few ambiguous elements as possible 
into the triangulation we are motivated to exclude that a v(A) is an interiour point of an 
edge or that it anyway does not contribute. The latter can be achieved by assuming that 
the edges Sj(A) have co-linear tangents at v in this case. 

Now we are left with those v(A) that are vertices of 7. Following the principle that our 
prescription should be canonical we must have that either each vertex of 7 is a basepoint 
of some A or none. Since the latter possibility is excluded by the principle of non- 
triviality we are now concerned with the issue of how many A's should have basepoint 
in each v G ^(7). A natural answer to this question is that there should be as many 
such A's as pairs of edges incident at v because otherwise we would single out one pair 
to another. However, we still need to fulfill the requirement that the A's must come from 
a triangulation. Both observations motivate to define a whole family of triangulations 
adapted to 7 and to average over them. 

Finally, we must fix in a diffeomorphism covariant way how to attach the arcs 012(A) 
to 7. Notice that since A is a part of a triangulation with v(A) a vertex of 7 and with 
Sj(A) segments of edges of 7 incident at v, it is possible that the endpoints of a 12 (A) are 
actually basepoints of of other triangles A'. This we either must avoid by choice (which 
is possible) or we must impose that the tangents of Sj(A) and 012(A) are co-linear at the 
endpoints of 012(A). As we will see, only the latter possibility leads to an anomaly-free 
theory. This furnishes our preliminary investigation of how to choose T(^y). 

We will now prescribe T(7). The prescription is simpler but very similar to the three- 
dimensional case. 

Fix a vertex v of 7 and let n denote its valence. We can label the edges of 7 incident at 
v in such a way that 

1) the pairs (ei, e 2 ), (e 2 , e 3 ), .., (e n _i, e n ) are right oriented and possibly also (e n , ei) is 
right oriented according to definition (|4.1|) and 

2) as one encircles v counter-clockwise, one does not cross any other edge after one crosses 
ei and before one crosses e$+i where e n+ i = e±. We are going to construct a triangle A 
associated with each such right oriented pair which we will call (ei,e2) from now on. 
We do not, in contrast to the 3+1 theory, construct a triangle associated with each pair 
because then oi 2 in two dimensions would intersect not only si,s 2 but also other edges 
of the graph which we must avoid in order to have an anomaly-free theory as we will see. 
Moreover, in two dimensions the way we ordered the edges incident at v is very natural 
and not available in three dimensions. 

Finally, let E(v) equal n if (e n , e\) is right oriented, otherwise let it equal n — 1. In 
particular, for n = 2 we must have E(v) = 1. 

We choose now S;(A) to be any segment of which does not include the other endpoint 
of Ci different from v and which starts at v. Furthermore, connect the endpoints of si(A) 
with the endpoint of 52(A) by an arc 012(A) with the special property that the tangent 
of oi 2 (A) is 

1) parallel to the tangent of Si(A) at the end-point of Si(A) and 

2) anti-parallel to the tangent of s 2 (A) at the end-point of s 2 (A). 
Two remarks are in order : 

a) Notice that we do not have to worry about any other edge of 7 intersecting 012(A) 
because in two dimensions the topology of the routing of 012 through the edges of 7 is very 
simple : there is no way that a 12 can intersect any other edges of 7 other than Si, s 2 given 
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the labelling of made above. This is in contrast to the three-dimensional case where 
the topology of the routing was extremely complicated to prescribe in a diffeomorphism 
covariant way. 

b) In contrast to the three-dimensional case we here prescribed the C 1 properties of the 
edges si, s 2 , a\ 2 at their intersection points. The reason for this will become evident only 
later when we prove anomaly-freeness. We will see that the C 1 property of the intersection 
is crucial. 

Whenever (e n , ei) is a right oriented pair the n triangles saturate v. Otherwise there 
are only n — 1 triangles and they do not yet saturate v. We follow the approach proposed 
in in order to achieve saturation. Namely, we take each of the E(v) triangles and 
construct three more from it such that they altogether saturate v. Then we average over 
the E(v) triangulations based on using only one such quadrupel of triangles. The details 
are as follows : 

Let Sj(t), ai 2 (t) be a parametrization of Sj, a% 2 with t G [0, 1]. Let 

Si(t) :=v- (Si(t) -v) = 2v- Si(t), 

a l2 (t) := 2v - a 12 (t), 

a- 21 (t) :=s 2 (l)+t( Sl (l)-S2(l)), 

a 2l (t) :=s 2 (l)+t(si(l)-s 2 (l)). 

Then it is easy to see that (si,s 2 ), (s 2 ,Si), (s 2 ,si) are right oriented pairs and that the 
four triangles A 12 , Ajj, A 21 , A 2 j based on these triples of edges saturate v (use ai 2 (0) = 
Si(l),ai2(l) = s 2 (l) to see this). 

Let now Si(v) denote the region in E filled by these four triangles based on a pair 
of edges (ei,e i+ i) incident at v. Also denote by Aj(v) the original triangle defined by 
s\, s 2 , d\ 2 for that pair from which we constructed the remaining three triangles as above. 
Let S(v) := uf^i'Si(v) be the union of these regions given by all the E(v ) pairs and let 
Si (v) = S(v) — Si(v). We will choose all the triangles so small that the S(v ) are mutually 
disjoint. Finally, let S = U„ e y( 7 )5'(u) and S = E — S. Then we can trivially decompose 
any integral over S as follows 




In the last line we have exploited that for smooth integrands and small triangles the 
integral over each of the four triangles constructed is the same up to higher order in the 
parameter 8 introduced before equation ( fL5|) . It is clear that the term in the first square 
bracket of the last line in (|4.9| ) is a sum of integrals over regions of E each of which does 
not contain vertices of 7. 

We are now ready to specify the family of triangulations T{^) of E which by Q4.5Q 
can actually be reduced to a family of triangulations of S, Si(v), Aj(t>) for v G V(j),i = 
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l,..,E(v) : 

1) Triangulate Aj(t>) by A;(-u) 

2) Triangulate S and Si(v) arbitrarily subject to the condition that no basepoint of a 
triangle should lie on an edge of 7 or that all tangents at an intersection with an edge of 
7 are co-linear 

3) The triangles Aj(t>) collapse to v as T — > 00 in such a way that all graphs 7 U Aj(t>) 
are diffeomorphic as T — > 00. In fact as long as we keep the prescription of how to choose 
Si (A), 012(A) specified above, all the graphs 7U A are related by an analyticity preserving 
smooth diffeomorphism no matter how "large" A. Namely, such diffeomorphisms can leave 
the image of 7 invariant while putting ayi in any diffeomorphic shape. 

Notice that now we have a well-defined prescription for the continuum limit because by 
construction the triangles that triangulate S, Si(v) do not contribute to the operator (f4.8|). 
The fact that the number of triangles that have their basepoint in vertices of the graph 
(which are the only ones that contribute) stays constant (namely E(v)) indicates that the 
continuum operator will be densely defined. 

4.3 Continuum Limit 

Let us summarize : having specified the triangulation we have triangles A(7, T) associated 
with the graph, more precisely E(v) for each vertex v of 7, the index T indicating that 
the continuum limit has not been taken yet. Then the regulated operator (|4.8|) becomes 

H T (N)f 7 := # 7i t(A0/ 7 := i £ N ( V )(^~J E ^ x 

11 veV(y) ^V 3 ) v(A),v{A')=v 

X ^(^(AoVCAjI^A)' V^]^(A)[^(A)> ( 4 - 10 ) 

where we have dropped the dependence of the A on 7, T. Now, since as T — ► 00 all 
holonomies approach unity, the limit T — > 00 does not have any meaning on the Hilbert 
space Tt = L 2 (A/Q, dfio). Indeed, on smooth connections we would get zero while on 
distributional connections the limit does not exist. Thus, the limit T — > 00 must be 
understood in another way. Indeed, recall that we wanted to impose the Hamiltonian 
constraint actually on diffeomorphism invariant distributions ^ G Now, the operator 
Ht(N) defines for each T an operator (Ht(N))' on $' by the equation 

[(H T (N)Y^]((j)) := y(H T (N)(j)) W G (f) E $ (4.11) 

because Ht(N) has domain and range in $ which is dense in 7i, for each T. Now, if ^ is 
diffeomorphism invariant then 

[(H(N))>*)(fJ := lim [{H T (N))'*]{fJ = *(H^ To (N)fJ (4.12) 

1 — >oo 

for each function / 7 cylindrical with respect to a graph 7 and for each 7. In other words, 
the number \E f (if 7j r(iV)/' 7 ) does not change under variation of T which by prescription 
corresponded to a diffeomorphism and so on diffeomorphism invariant states we may 
evaluate it on any finite value T and the T — > 00 limit is trivial. It follows that on 
diffeomorphism invariant states the continuum limit is already taken for (Ht(N))'. 
In fact, it is easy to see that this result can be extended to any product 
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[Ht{Ni)Ht{N 2 ) ..HxiNn)]' because the triangles attached have, at each level, an unam- 
biguously defined diffeomorphism covariant location. This observation is needed in order 



to give sense to commutator computations P, |T2|. 

In the sequel we will drop the index T and understand that when finally evaluating ev- 
erything on diffeomorphism invariant distributions the value of T is irrelevant. 



5 Consistency 

There are two kinds of consistencies to be discussed : 

The first is the cylindrical consistency, that is, we have obtained a family of operators 
(if 7 (iV)) 7 which should be projections to cylindrical subspaces of a "mother" H(N). That 
such a H(N) exists has to be proved. 

The second is that we need to make sure that H(N) does not suffer from quantum 
anomalies. 

5.1 Cylindrical Consistency 

In proving that a family of operators (Oj, Di) ie i on a Hilbert space 7i, where Di is the 
domain of Oj and where I is some partially ordered index set I with ordering relation <, 
is cylindrically consistent we need to reveal that whenever i < j that Oj is an extension 
of Oj, that is 

1) The domain of Oj is contained in that of Oj, D L C Dj and 

2) The restriction of Oj to Di coincides with Oj, {Oj)\i = Oi. 

Let us check that this is the case for our operator family. Recall that a spin-network 
state depends on all of its edges non-trivially in the sense that all edges carry spin j > 0. 
The space <3? 7 is the set of finite linear combinations of spin-network states which depend 
on the graph 7. Now, while the set of graphs can be partially ordered by the inclusion 
relation, the set of cylindrical functions cannot because a function which is defined on 
a smaller graph is defined also on any bigger graph that properly contains it, however, 
the additional edges in that graph automatically carry spin zero and so the cylindrical 
subspaces cannot be compared. Another way of saying this is that given a cylindrical 
function / we can uniquely decompose it as / = XL / 7 , / 7 e $ 7 and on / 7 we have 
unambiguously H(N)f^ = H^(N)f^. We cannot write H(N)f~, = H^(N)f^ with 7 C 7' 
because there is a condition on the spins of the edges of 7' involved when applying Hy 
which is not satisfied for / 7 . In other words, $ 7 D $y = if 7 7^ 7' 

We conclude that the family (if 7 (iV)) is trivially cylindrically consistently defined and 
therefore defines a linear operator on all of TC. 

5.2 Anomaly- freeness 

Recall that the classical Dirac algebra is given by 

{H(M), H(N)} = d 2 x(M, aN - MN, a )q ab V b 

where V a is the vector constraint. That is, the Poisson bracket between two Hamiltonian 
constraints evaluated on the constraint surface defined by the diffeomorphism constraint 
vanishes. 
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In the quantum theory one would therefore like to verify that naively 
[H(M), H(N)]f = for any state / that satisfies V a f = 0. Several subtleties arise : 

1) The solutions V a f are in general no elements of the Hilbert space but generalized 
eigenvectors (distributions). Indeed, in this context the solutions of the diffeomor- 
phism constraint are not elements of 7i but of $' where we have the proper inclusion 
$ C H. C Thus, since H(N) is defined only on $, the only operator that is de- 
fined on $' is the dual (H(N))' via the pairing V(H(N)<p) = [(H(N))'^]((p). 

2) Observe that the operator (H(N))' was not defined on every distribution but actu- 
ally only on those that are solutions of the diffeomorphism constraint. Now even 
if ^ is diffeomorphism invariant, that is, ^/(Lf^ip)^) := ^(f^t^) = ^(/-y), then 
(H(N))'ty is not any longer as one can easily check. Thus we cannot verify that 
[(H(M))', (H(N))']ty = 0, this equation is simply not defined. However, what is 
well-defined is ([H(M), H(N)])'fy = and this is what we are going to verify. In- 
deed, there is no hope to make sense out of [(H(M))', (H(N))']^ since not even 
classically H(M) is diffeomorphism covariant. On the other hand, one could pro- 
ceed as in |12]| and define H'(M)H'(N) := (H(N)H(N))' which makes sense again 



on diffeomorphism invariant states. 
3) One might be even more ambitious and ask that 



([H(M), H(N)})' = ( / d 2 x(M >a N - MN a )q" b V b )' (5.1) 



s 

that is, the Dirac algebra is faithfully implemented in the quantum theory. However, 
there are several issues that prevent us from doing so. First of all, the generator 
of diffeomorphisms, V a , does not have a quantum analogue, the diffeomorphism 
group does not act strongly continuously on H. So the only thing that we can 
hope to obtain is something like O f [U((p) — 1] for the right hand side of ( |4.12|) for 



some (p £ Diff(S) and some dual operator O' (notice that U(<p)' = U(cp ) can 
be extended to all of $'). Secondly, the situation is even worse for q ab . Thirdly, 
since, as we said, ([H(M), H(N)})' is only well-defined on diffeomorphism invariant 
distributions, then either the dual of the commutator vanishes or it does not. In the 
latter case there is an anomaly even in the sense of U{4>) — 1. In the former case we 
get just zero but then we can trivially make an equality of the form 

([H(M),H(N)]y = d'[U(<p)-l] 

for any O' that we like. It then remains to ask whether one can somehow make sense 
out of an operator corresponding to the combination J s d 2 x(M <a N — MN !a )q ab Vb and 
that is actually the case : We will not prove this assertion here but refer the reader 
to WM which treats the 3+1 case but from which it is obvious that the result can 



be extended to the 2+1 case. 
Summarizing, we will check that QiT(M), H(N)]y^ = on diffeomorphism invariant 



states. The key element of the proof is the following : as is obvious from ( 4.10 ), if / 7 is 
a function cylindrical with respect to a graph, then H(N)f 1 is a linear combination of 
functions each of which depends on graphs with new vertices not contained in 7. More 
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precisely, if v G V{pj) then for each triangle A based at v there is a term ^h:H A f n an d 
this function is a linear combination of functions /' each of which depends on a graph 7' 
contained in the following list : 

7 U A, (7 U A) - si(A), (7 U A) - s 2 (A), (7 U A) - (s a (A) U s 2 (A)). Whether they appear 
depends on the spins of the graph 7. In any case these functions /' depend on two more 
vertices v±, v-i coming from the endpoints of the arc 012(A). They may not depend on the 
original vertex v if that vertex was two-valent with spins of the edges corresponding to 
Si (A) being j = 1/2 for both i = 1,2. In that case [h~} A r v \\, V v ]f = because neither /' 
nor h~r A r v \->f depend on graphs with more than one edge incident at v. 
The point is now that [h~^ A ( Vi ^, V Vl ]f = [^(A(« 2 ))' ^a]/' = 0- The reason for this is that 
the vertices v±,V2 in the graphs on which /' and hj5 A ( v \\f depend does not have edges 
with linearly independent tangents incident at it so that the volume operator annihilates 
these functions. 

Let us now write ( f4.1U| ) in the form 
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Hy(N) = T2 E N(v)H 7tV 

11 vEV(-y) 

H-y,v = p/..\2 E H-r,v,AA' 

rjy - U ' v(A),v(A')=v 

Hy,vAA' = ^ 3 ^to( h o^(&)h Sk (A)[h^ A ),yfyh^^ (5.2) 

The function H lv f 1 now can be written as a linear combination of functions f , each of 
which depends on a graph 7' which is a proper subgraph of the graph 7(f) := 7LUa)=u(A')=u 
[A U A'] and we will mean by H-y^y the operator that reduces to Hy tV > on fy for each 
v' E V(j') and is zero if v' V(y). 
With this preparation we evaluate 

[H(M),H(N)]f 7 = Yl [N V H{M) - M V H{N)\H^ 

veV(-y) 

= E E [N V M V > — M v N v t]Hy v -) )V i H ljV j 1 

i)6V(7) v'eV(y(v)) 

= ^ [N V M V > — M v N v i]Hy( v ^ v tH~f )V f~f 
v,v'ev(j) 

= 2 E [iV»M«' - M^iVj/J^^yEy^ - Ey^/)^^^/]^ . (5.3) 

v,v'GV(-f) 

Here we have used our notation to write the commutator as a double sum in v G V(j),v' G 
V(l( v )) i n the second line, then in the third line we have used the important fact that 
the constraint does not act at the new vertices that it creates so that the sum over 
v' G V(j(v)) collapses to a sum over the original v' G V(j) and in the last step we have 
used the antisymmetry in the lapse functions to write the product of operators as their 
antisymmetrized sum of products. Clearly the term with v = v' vanishes trivially. If 
v ytz v ' then, since H 1)V manipulates the graph only in a small neighbourhood of v, we can 
commute the two operators in the last line of ( |5.3|) to write both with the vertex v to the 
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right hand side as 

[H(M),H(N)]f 7 — ~ J] [^M„,-M„iV„,][^ 7My ^-F 7y ^ 7( „ 0) „]/ T . (5.4) 

»,t)'6V(7) 



Now by inspection of ( |5.2j ) we see that the last square bracket is a linear combination 
of functions of the type f — f where /, /' are related by an analyticity preserving dif- 
feomorphism by construction of the triangulation which relates different choices for the 
loop attachment by such a diffeomorphism (this point is explained in more detail in ||). 
Thus when evaluating ( |5.4j ) on a diffeomorphism invariant state we can remove those dif- 
feomorphisms and obtain just zero. 
This suffices to show ([H (M) , H (N)])' = 0. 

Notice that it was essential in the argument that the additional vertices created by 
H{N) when acting on / 7 do not contribute as we showed. If that was not the case the 
commutator would not vanish on diffeomorphism invariant distributions which is why we 
attached the loop in such a particular, C 1 , way. 

6 Solving the theory 

This section is divided into two parts : In the first part we will describe the complete space 
of solutions to both the diffeomorphism and Hamiltonian constraint. In the second part 
this solution space is shown to contain the solutions to the curvature constraint which 



can be formulated in our language as well JT3[. One can equip the solution space with at 



least two very natural inner products. One of them is the inner product appropriate for 
the curvature constraint, the other one arises from direct construction of the solutions in 
the first part this section. Neither of these inner products give all solutions a finite norm. 

6.1 Complete set of solutions to all constraints 

Let be given a spin- network state T 7 j s and let 

i T ,J,c} ■= &(<PKfc V G Diff(S)} 

be its orbit under the diffeomorphism group of analyticity preserving smooth diffeomor- 
phisms. We define a diffeomorphism invariant distribution on $ by 

That this is a continuous linear functional on $ follows from the fact that the spin-network 



states form an orthonormal basis by the argument given in |TI|. Therefore It is also clear 

that every diffeomorphism invariant state is a linear combination of such [T j -] 's so that 

by this procedure we can claim to have found the general solution \l/ to the diffeomorphism 

constraint W(U((p)f) = *(/) VV G Diff(E). We will call this space $' Di// 

Given any / G $ we can uniquely decompose it as / = J2i fiTi where fj are some 

constants and Tj are spin-network states. We then define [/] := rjoifff '■= 

Notice that one cannot define [/] as the sum of all states which are in its orbit under 

diffeomorphisms since spin-network states defined on different graphs have uncountably 
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infinite sets of diffeomorphisms that move one graph but not another. We have been here 
imprecise with the issue of graph symmetries which alter the above formulae somewhat. 



See |T2[ for more details. 

This construction can be used to define an inner product on ^'oiff by 

<[f],[g] > mff := [f](g) 

which is clearly a positive definite sesquilinear form and equips $£>i// w ith the structure 
of a pre-Hilbert space. 

Remark : 

It has been shown in |]] that if there are only strongly diffeomorphism invariant ob- 
servables in the theory then those observables define a superselection rule, namely, they 
cannot map between spin-network states based on graphs which are in different diffeo- 
morphism equivalence classes. As a result, the group average could be defined differently 
in every sector, that is, the inner product in every sector can be chosen individually which 
amounts to the ambiguity that the particular way of averaging given by [/] is not selected 
by physical principles, meaning that for every diffeomorphism equivalence class of graphs 
there could be a different constant that multiplies [TL7J in [/]. 

However, as there are clearly weakly diffeomorphism invariant observables which, together 
with the Hamiltonian constraint map between those sectors, there is no superselection rule 
and the way we have averaged is selected by the requirement that averaged spin-network 
states remain orthonormal [O. 

We now wish to employ this result to find the general solution to all constraints. To 
that end, consider the set 

R:={H(N)(j), N eS, 0G$}, 

the range of the Hamiltonian constraint on $ where S denotes the usual Schwartz space 
of test functions of rapid decrease. Consider its orthogonal complement in $ denoted 
S := R L C $. Finally, consider the set {[s], s G S}. Then it is easy to see that every 
solution to all constraints is a linear combination of elements of this set and we will call 
the resulting span $' phys . Namely, let s = J2i s iTi G S then by definition J2i $ifi — for 
any / = £j fjTj G R. Thus [*](/) = £i sj[Ti](f) = Ej hfi = OV/ G R. ' 

A geometrical construction of the space S was given for the three-dimensional the- 
ory in ||. Here we could proceed similarly. However, since this is only a model we 
restrict ourselves to showing that the space <&' phys is uncountably infinite dimensional. 
Namely, a particular simple class of vectors in S consists of those elements of $ which 
are linear combinations of spin-network states whose underlying graph is not of the form 
7 Ua 12 (A), [ 7 Ua 12 (A)] - Sl (A), [ 7 Ua 12 (A)] -s 2 (A), [ 7 Ua 12 (A)] - [s a (A)Us a (A)] for any 
A = A( 7 ), v (A) G V( 7 ) and 7 is a graph underlying the same restriction but is otherwise 
arbitrary. This particular class of solutions has the property that all of the resulting [s] 
are normalizable with respect to < ., . >Diff while genuine elements of &' phys will not 
be normalizable with respect to the kinematical inner product on ^'mff- O n the other 
hand, since the 2+1 Hamiltonian constraint really resembles the 3+1 Euclidean Hamil- 
tonian constraint it follows from the redults on the kernel of the Euclidean Hamiltonian 
constraint given in 0, |12| that every solution is a (possibly infinite) linear combination of 
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basic solutions each of which is in fact normalizable with respect to < ., . >Diff- 
We will see that the solutions to the curvature constraint are not normalizable with respect 
to < ., . >Diff and one needs to define another appropriate inner product < ., . > cn ™ on 
the subset of Q' phys corresponding to the solutions of the curvature constraint. However, 
it will turn out that the natural inner product < ., . > p h ys for our Hamiltonian constraint 
as suggested by [12] is such that curvature constraint solutions are still not normalizable 
and so < ., . > cur v and < ., . > p h ys define genuinely non-isometric Hilbert spaces. We will 
turn to that issue in the next subsection. 



6.2 Comparison with the Topological Quantum Field Theory 



As shown in |13[], in our language a solution to the curvature constraint = in the 
quantum theory is a distribution $j G $' given by \I/ / := 5 fJj0 (F)f for any / G $. Here 
S fi0 (F) is a 5 distribution with respect to the inner product on 7i which has support on 
the space of flat connections modulo gauge transformation M.. More precisely, we have 
the following : Any function on Ai is a gauge invariant function which depends on the 
connection only through the holonomies along (representants of) the independent gener- 
ators ai,.., a n of the fundamental group 7Ti(E), that is, f(A ) = f n (h ai (A ), .., h an (A )) 
for A G M.. The measure z/ on M. for gauge group G is defined by 



duo(A Q )f(A Q ) := / d/2 H (g l )..dno(g n )f n (g l} ..,g n ) 

M JG n 

where /j,h denotes the Haar measure on G. Then the delta distribution for flat connections 
is given by 



S W (F(A)):= dp (A )5 IM) (A ,A) (6.1) 

J M 

where 



6, (A , A) = £ T l3 M)T l3 ,{A Q ) (6.2) 

and the sum runs over all possible spin-network states. It is possible to arrive at (|6.1|) 
from first principles by following the group average proposal [13|. 



It is also possible to write ( |6.1|) as a linear combination of distributions in ■ To that 
end, denote by / the label of a spin-network state and define Tj/j := [Ti]. Notice that the 
integral kj := J M dv^A^T^Ao) =: k[r\ is diffeomorphism invariant and thus only depends 
only on [/]. Then we may write 



5 fl0 (F(A))=Y,k [ i] T U]( A )- (6-3) 

m 

It is easy to see [0 that (|6.3|) is a distribution on ^'aff and certainly it is a distribution 
on However, ( |6.3| ) is not normalizable with respect to < ., . >Diff '■ 
To see this we use (|6~3| ) to notice that we can write 5(F(A)) = r] Di fff F where fp '■= 
J2\i] k[i]Ti ({i])(A) and -ZoQ-TI) e [I] * s an arbitrary choice. Thus by definition of the inner 
product between diffeomorphism invariant distributions we find 

\K(F)\\ 2 Diff = ( V DifffF)(fF) = E \ k m\ 2 (6-4) 

[I] 
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where the sum is over diffeomorphism equivalence classes of spin-network labels. But 
quantity (|6.4j ) is just plainly infinite : 

Namely, it follows from the definition of kj = k\n that k\n = Tj(A = 0) whenever the 
graph underlying / is contractable. There are an at least countably infinite number of 
contractable, mutually non-diffeomorphic, non-trivial graphs 7„, n = 1,2,.. in any S. 
An example is given by choosing 7 n to be an n-link, that is, a union of n mutually non- 
intersecting loops ati, .., a n homeomorphic to a circle each of which is homotopically trivial 
(contractable). Choose I n such that T In (A) = Y\H = iT ak (A) where T a (A) := ti(h a (A)) is 
the Wilson-Loop function and h a (A) denotes the holonomy of A along the loop a. Using 
the basic integral f SU ( 2 ) d\iH{g)9AB9cD = \^ac^bd it is easy to see that Tj n provide an 
orthonormal system of spin- network states. But Tj n (A = 0) = 2 n and so ( |6.4|) contains 
the meaningless sum Y^=i 2 n - 

We must check whether or not ^/ is also a solution to the constraint H(N) (it obvi- 
ously is diffeomorphism in variant). To that end we must compute 

9f(H(N)fJ = [ d^5, (F(A))(fH,(N)f,)(A) 

= I du (A )(fH 1 (N)f 1 )(A ) = (6.5) 

J M 

because either H y (N)fy is identically zero or it is a linear combination of the vectors 
(recall Q) 

H 7jVj A,A'f-r = ~ 2e i - J e fc/ tr(/i a . j .( A ')T m )tr(r m /i Sfc ( A )[/i^ 1 (A) , \[^v]h Sl (/±)[h~^ A - ) , \ &v))f 7 

which therefore are proportional to the matrix elements of [h ai2 — h~* ] for a contractable 
loop a which vanishes on A e M.. Here we have used the su(2) Fierz identity tr(rjA)tr(rj-B) = 
tr(A)tr(S)/4 - tr(AB)/2 together with e ij ti(h aij ) = tv(h ai2 ) - tr(h~* 2 ) = 0, a particular 
property of SU{2) (we did not need to use this, the result holds for general G). 
Thus, any solution to the curvature constraint is a solution to the Hamiltonian constraint. 
However, as we have demonstrated, there are an infinite number of more solutions to 
the Hamiltonian constraint, in particular those which are normalizable with respect to 
< ., . >Diff an d no solution to the curvature constraint has this property. Notice that the 
inner product on the space of solutions to the curvature constraint comes from a group 
averaging map, it is just given by (|L3| ) 



< V f ,V g > C urv-= Vf(g) = 




It is now tempting to view this result as the restriction to the special solutions of the 
curvature constraint of a more general inner product appropriate for the Hamiltonian 
constraint. 

There seems to be an unsurmountable obstacle : the Hamiltonian constraint is not a 
self-adjoint operator on TC and so group averaging as defined in || cannot be employed. 
Moreover, group-averaging really means to exponentiate the Hamiltonian constraint and 
that in turn implies that we know the motions it generates and thus we would have to 
completely solve the theory. Thus, it seems that we cannot define a map 
^'phys'i f ~^ Vf- However, in the case that we have self-adjoint constraint operator, 
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the group-average algorithm is nothing else than a sophisticated way to construct the 
projector onto the distributional kernel of the constraint operator (this is explained in 
more detail in |I2||). We are therefore lead to define the map rj, in the case that we do not 



have a self-adjoint constraint operator, as a certain (generalized) projector on the kernel 



of the constraint operator. As in [13] we split the problem into two parts and proceed as 
follows : 

Given / G $ we have a group averaging map i] Di ff : $ — > &Diff defined by rjDiff(f) '■= [f] 
and an inner product defined by < [/], [g] >mff-= L/](sO : =< [f\,9 >■= fj/gd/M)[f]9- We 
define now $Ham '■= &Diff an d would like to define a map T]Ham '■ &Ham — ► ^Ham- The 
space Q' Ham coincides with <&' phys when viewed as a space of distributions on $ via the 
map r] := rjHam ° VDiff- It remains to construct r\uam- 

As we have seen, the elements [s] G & p h ys , s G S span Q' phys . Moreover, by explicit 
construction (given for the 3+1 theory in H) we can orthonormalize them with respect 
to < ., . >mff thus exploiting that in our case all these [s] are normalizable with respect 
to < ., . >Diff. We obtain particular elements G ®Ham ^ ^Ham with the property that 
< if>p,if>„ >Diff= &w We are now ready to define the projector t] Ham : given if) G § H am 
define 

VHam^ ■= VVVVW : = < ^ ^ >Diff ' ( 6 ' 6 ) 

Notice that even if not all of the [s] would be normalizable with respect to < . , . >mff 
then one could still take (|6.6|) as the group average map, just the elements ip^ now form a 
basis in the generalized sense that they are mutually orthogonal in the sense of generalized 
eigenvectors (similar to usual momentum generalized eigenfunctions of ordinary quantum 
mechanics which are not really orhtonormal in the Hilbert space sense but only orthogonal 
in the sense of 5 distributions). 

The fact that the Vv are normalizable with respect to < .,. >Diff displays r\uam as a 
projector on a genuine subspace of Hmff- 

Observe the dual role of the which we can view both as elements of $# am and as 
elements of &Ham = &Diff ^^Diff- in particular, notice the peculiar identity rj Ham ifj^ = 

We now simply define an inner product on the elements rjHamip by 

< VHam4>,VHam4>' >Ham- = fa/Tom == ^2 < VV> ^ >Diff < VV> >Diff (6.7) 



for each if), if)' G §Ham- Expression (|6.7|) is clearly a positive semi-definite sesquilinear 
form with the property that the 77^ remain orthonormal. It is also independent of the 
orthonormal system if)^ (the label /1 is "nicely split" into a discrete piece and a continuous 
piece and ip^s are orthonormal with respect to both pieces in the sense of Kronecker S's, 
see 



12 1 for details). 



We now combine the two group average maps to obtain 

V ■= Vphys ■= VHam ° VDiff ■ $ := %hys -> % hys ■= $ff am , / -> VHam[f] = 

ft 

(6.s; 

and the physical inner product for the elements rj p hy S f becomes 

^ Vphys f ', 1]phys f ^phys *^ T) H am[f\ , 1] H am [9] > Ham 
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J2 < ^ if] >Diff < id] >Diff 



= J2Mf)Ma) = (v P hy S f)(g) (6.9) 

where in the second before the last equality ip^ is viewed as an element of 

Notice that with this definition H p h ys C Hniff- That this makes sense is shown in fT2| . 



In other words, infinite linear combinations of elements of are allowed but only with 
suitably converging coefficients. 

As we have already shown that ^/ ^ Hdiff it follows that ^ 'Hphys-, no solution to 
the curvature constraint is normalizable with respect to < .,. >Ham- Since what really 
determines the physical inner product is the Hamiltonian constraint, for instance via the 
group average approach, this result was expected given the totally different algebraic 
structure of the two sets of constraints. In particular, the scalar product < ., . > cur „ is 
rather unnatural from the point of view of the Hamiltonian constraint. 
The reverse question, whether 1 1 . 1 1 jjam normalizable elements of $# am have finite norm 
with respect to ||.|| cun ; cannot even be asked in general because a general element of &' Ham 
cannot be written as tyf, / G $. 

We conclude that the sectors of the theory described by either of the inner products 
< •,• >Curv and < .,. >ph ys are mutually singular (that is, the underlying measures of 
the scalar products are singular). On the other hand, as far as the space of solutions 
to the constraint is concerned we find that all solutions to the curvatur constraint are 
annihilated by the Hamiltonian constraint. Moreover, if we choose < .,. >curv as the 
inner product then we find complete agreement with the results in [ffl, ^ although our set 
of constraints and our quantization approach was totally different from the outset. Since 
we copied step by step the quantization procedure of || to maximal extent, we conclude 
that this procedure does lead to the correct answer in the present model which is a small 
but non-trivial check whether the proposal of fU is reliable or not. 



7 Conclusions 

The aim of the present paper was to check whether the method of quantizing 3+1 general 
relativity by the method proposed in [Q] is reliable in the sense that when that procedure 
is applied to well-known models we get the known the results. When applied to 2+1 
Euclidean gravity we find complete agreement thus giving faith in those methods, the 
more, as 2+1 Euclidean gravity is maximally similar to 3+1 Lorentzian gravity as far 
as the algebraic structure of the constraints and the gauge group are concerned (at least 
when we consider the Hamiltonian rather than the curvature constraint). 

On the other hand, the quantum theory as obtained by our approach has a much bigger 
space of solutions to all constraints than the space as obtained by traditional approaches 
while the latter is properly included in our solution space. A natural question that arises 
is then what to do with those extra solutions and how to interprete them. In particular, 
there seems to be a clash between the number of classical and quantum degrees of freedom. 

Now, a hint of how to interprete these solutions is that many of them are of the form 
[s], s 6 S and so s is a cylindrical function. Therefore the volume operator V(B) van- 
ishes on s for almost every B. This suggests that [s] is a spurious solution because if 
we want the classical theory defined by curvature and Hamiltonian and diffeomorphism 
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constraint to be equivalent (recall that we had to impose det(g) > classically) then we 
must really ask that the volume operator V(B) is strictly positive for any B, before tak- 
ing the diffeomorphism constraint into account. We conclude that no cylindrical s should 
give rise to an element of $' phys via s — > [s] (which can be achieved by superposition of 
an infinite number of the [s] or by considering infinite graphs) which would presumably 
remove the clsh between numbers of degrees of freedom alluded to above. This latter 
observation gives rise to the speculation that also many of the solutions found in || for 
the 3+1 theory should be spurious because in the classical theory we have to impose the 
anholonomic condition det(g) > as well. On the other hand it may be desirable to allow 
for degenerate solutions at the quantum level because by passing through singularities 
of the metric one can describe changes in the topology of the hypersurface S. Therefore 
one may expect that some of the solutions actually carry topological information and, 
moreover, that although we have started from a fixed topology in the classical theory 
we end up describing all topologies at the quantum levelQ. The complete answer to this 
puzzle is left to future investigations. 
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A Spectral analysis of the two-dimensional volume 
operator 

Since the gauge group is still SU (2) we may copy the results from || to compute the full 
spectrum of the two-dimensional volume operator. In particular it follows immediately 
that this operator is essentially self-adjoint, positive semi-definite and that its spectrum 
is entirely discrete. This holds on either gauge invariant or non-gauge invariant functions. 
In this appendix we restrict ourselves to the part of the spectrum coming from graphs with 
vertices of valence not larger than three, that is, we display the eigenvalues of the operator 



V v of (|3.13|) restricted to vertices v of valence n = 2, 3 on gauge invariant functions. Indeed, 
as the volume operator cannot change the graph or the colouring of the edges of the graph 
with spin quantum numbers of a spin-network state it follows that it can change at most 
its vertex contractors. However, given the spins of the edges incident at v, the space of 
vertex contractors is one- dimensional for n = 2, 3 by elementary Clebsh-Gordon theory. 
Therefore spin-network states all of whose vertices have at most valence three must be 
eigenvectors of the volume operator (in any dimension). Notice also that all the V v for 
different v's are mutually commuting. In three dimensions these spin-network states are 
in the kernel of the volume operator, in two dimensions none of them is annihilated as we 



2 This speculation on the conceivably topological meaning of the solutions is due to Abhay Ashtekar. 
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will show (as long as the tangents of the edges at v span a plane). 

For n = 3 there are only two generic non-trivial situations : Either (Case A) no two 
of e%, e-ii 63 have co-linear tangents at v or (Case B) two of them, say e%, e 2 have co-linear 
tangents at v but not ei, or e 2 , e%. Here e%, e 2 , 63 are the three edges incident at v which 
are coloured with spins 31,32,33 e {ji +32,31+32-1,-, \ji ~ h\} respectively. We can get 
the eigenvalue for the case n = 2 by taking the result for n = 3 and setting for instance 
3s = 0,ji =32 = j + 0. 

In the calculations that follow we will use the following notation : 
Qv ■= [^V v ] 2 = KEl 

K := \ £ sgn(e /? e J )X} J where Xjj := e ijk XjX$ 

1</,J<3 

X] := Xi r , Xj := (X)), X u = X]X), Aj := X n (A.l) 

and it is implied that J, J,K e {1, 2, 3} are mutually different so that [X], Xj] = 0. As 
the notation suggests, A/ is the Laplacian on SU{2) with spectrum + 1), 2j > 
integral. Notice that X % u = —Xjj so that 

K= £ sgn(e 7 , ej )^j . 

1</<J<3 

As in the main text we will use generators of su(2) with structure constants +€ijk which 
implies that [X}, Xj] = —eijkXf and so eijkX}Xj = —Xf (the minus sign comes from the 
right rather then left invariance) . 

There are some identities among these quantities that we are going to use. The first one 
is the familiar spin recoupling identity 

2Xjj = [Xj + Xj] 2 -A I -Aj = A K -Aj- Aj (A.2) 

where in the second equality we have used the fact that Xi + X 2 + X 3 = 0, that is, the 
total angular momentum operator vanishes of on gauge invariant functions. Then if / is 
gauge invariant 

[Xj + Xjff = — [Xi + Xj]X K f = -X K [Xj + Xj]f = [X K ] 2 f = A K f 

and of course the Aj commute with every X). The next identity is, using basic e^k 
arithmetic 

Xh = X}Xj(X}Xj - XjXj) = A/Aj - X}([Xj, X)] + X}Xj)Xj = A T Aj + X„ - X} 3 

(A.3) 

and by very similar arguments 

XjjXjk = —AjXjk + XjjXjk + €ijkX\X 3 K X k J . (A. 4) 

The last term in (|A.4 ) is essentially the basic operator from which the tree- dimensional 
volume operator is built and which vanishes in the three-valent case on gauge invariant 
functions. Indeed, replacing, say Xj = —Xj — Xk and using the su{2) algebra for the X\ 
we see that that term vanishes. 
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Remarkably, upon substituting for Xjj according to ( A.2 ) we find 



XijXjk = |[2(A 7 A J + AjA K + A K Aj) - (A? + A] + A^)] (A.5) 

which is independent of the choice of the pairs {U), (JK). 

We have now all tools available to finish the calculation. We will treat cases A, B 
separately. 

A) We may label edges without loss of generality such that sgn(ei,e 2 ) = sgn(e 2 ,e 3 ) = 
sgn(e 3 ,ei) = 1, that is, we cross ei,e 2 ,e 3 in this sequence as we encircle v counter- 
clockwise. Then E v = X 12 + X 23 + X 3i . We just need to use (|A.2| )- ([A~75| ) and to be 



careful with the order of J, J in Xu to find after tedious algebra that q v = [E v ] 2 is 
just given by 

q v = ^[2(A 1 A 2 + A 2 A3 + A 3 A 1 )-(A2 + A 2 ! + a2)] 

- ~(Ai + A 2 + A 3 ) . (A.6) 

Thus, the eigenvalue is obtained by replacing Aj by + 1). Expression (|2.6|) 

looks worrysome : is the eigenvalue going to be non-negative ? A moment of 
reflection reveals that it is even strictly positive unless ji = j 2 = j 3 = in which 
case it vanishes : It will be sufficient to show that the operator in the first line 
of ( |A.6|) has non-negative eigenvalue. We just need to remember that j 1; j 2 , j 3 are 
not arbitrary. We may assume without loss of generality that j 2 > ji such that 
3s e {ji + J2J1 + 32 - 1, ..,32 - 3\} ■ We have 

/(As) := 2(A X A 2 + A 2 A 3 + AgA^ - (A? + A 2 2 + A3) = 4A X A 2 - (A3 - A 1 - A 2 ) 2 

(A.7) 

which takes, in terms of eigenvalues, its lowest value at maximum value of the 
function IA3 — Ai — A 2 |. Given arbitrary j% < j 2 , since —A3 is a strictly increasing 
function of j$, we find that the extrema of that function are found for the extremal 
values j 3 = j 2 ± ji and are given by \2jij 2 \ and | — 2ji(j 2 + 1)| respectively. Then 
(E3 reveals that /( A a) > 4 Ji(j2 + l)(j2 - h) > because j 2 > j x . 



In case that we consider a two-valent vertex, we may just set A 3 = 0, Ai = A 2 = A 
and find the extremely simple result 

q v = —A . (A.8) 

B) We may, without loss of generality, label edges such that ei,e 2 have co-linear tan- 
gents at v (that is, sgn(ei,e 2 ) = 0) and such that sgn(ei,e 3 ) = sgn(e 3 ,e 2 ) = 1. 

Then E v = X 13 + X 32 . The same algebraic manipulations show that we get now for 

q v = [E v ] 2 the expression 

q v = [2(A X A 2 + A 2 A 3 + A3A1) - (Aj + A 2 2 + A3)] - A 3 (A.9) 

which is positive unless, of course, A3 = in which case it vanishes. 

We conclude that the two-dimensional volume operator has a much smaller kernel than 
the three-dimensional one, in particular, two and three-valent vertices, whether gauge 
invariant or not, do not contribute to the kernel. 
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